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ABSTRACT 

We present Subaru/IRCS J band data for Fomalhaut and a (re)reduction of 
archival 2004-2006 HST/ACS data first presented by Kalas et al. (2008). We 
confirm the existence of a candidate exoplanet, Fomalhaut b, in both the 2004 
and 2006 F606W data sets at a high signal-to-noise. Additionally, we confirm 
the detection at F814W and present a new detection in F435W. Fomalhaut b's 
space motion may be consistent with it being in an apsidally-aligned, non debris 
ring-crossing orbit, although new astrometry is required for firmer conclusions. 
We cannot confirm that Fomalhaut b exhibits 0.7-0.8 mag variability cited as 
evidence for planet accretion or a semi-transient dust cloud. The new, combined 
optical SED and IR upper limits confirm that emission identifying Fomalhaut 
b originates from starlight scattered by small dust, but this dust is most likely 
associated with a massive body. The Subaru and IRAC/4.5 upper limits 
imply M < 2 Mj, still consistent with the range of Fomalhaut b masses needed 
to sculpt the disk. Fomalhaut b is very plausibly "a planet identified from direct 
imaging" even if current images of it do not, strictly speaking, show thermal 
emission from a directly imaged planet. 
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Introduction 



In November 2008, two teams announced what were described as the first di rectly imaged 



extrasolar planets arou nd the nearby young stars HR 8799 and Fomalhaut ( jMarois et al. 



2008 



data ( iLafreniere et al 



Kalas et al.l 12008 ). Numerous studies recove r ed the three HR 8799 planets in archiva l 



20091 : iFukagawa et al.ll2009l : iMetchev et al.l 120091 : ICurrie et al.l 120121 ): 



recent ones identify a fo urth planet (HR 8799e, Marois et al. 2 010; Currie et al. 2011a) 
with a similar spectrum ( Skemer et al. 2012; Esposito et al. 20121 ). Detailed dynamical and 
atmospheric mo deling confirms that HR 8799 bcde are bona fide planets, not brown dw arfs 
{Mbcde < 9 Mj, iMarois et allboiOal : ICurrie et allboilat ISudol and Haghighipouil 120121 ). 



In contrast, refereed follow-u p studies ha ve yet to confirm Fomalhaut b, casting serious 



doubt on its status. As noted by iKalas et al.l . Fomalhaut b's spectrum looks little like that 
expected for a substellar object, as it i s detected in the optical but not in the inf ra red. Deeper 



IR data do not recover Fomalhaut b (jMarengo et al.ll2009l : iJanson et al.ll2012[ ). iKalas et al. 



describe Fomalhaut b's emission as being dominated by scattered dust emission surrounding 
a planet, with a therrnal com ponent at 0.8 fim and accretion-driven variability at 0.6 fim. 
However, iJanson et al.l (|2012l ) argue, based on their detection limits, that essentially 'none' of 
the optical emission originates from a planet, the 0.6 /im variability instead identifies a semi- 
transient dust cloud and thus Fomalhaut b is not a directly imaged planet . While Foraalhau t 
b seems to be a natural explanation for the well-confined debris ring (IKalas et al.ll2005l ). 
shepherding plane ts unassociated with Fomalhaut b may provide an alternate explanation 
teolev et al.ll2012h . 



These doubts about Fomalhaut b's status call for new data to further constrain its IR 
emission and a reanalysis of the HST data from which Fomalhaut b was originally identified. 
Thus, in this Letter, we p resent new Subaru/IRCS J band imaging of Fomalhaut and a 
rereduction of iKalas et al.l 's HST/ACS data. We use these data to constrain the near-IR 



emissi on from any planet corresponding to Fomalhaut b and revisit key results of iKalas et al. 



(120081): 1) Fomalhaut b's existence, 2) its spectrum, 3) 0.6 variability, and 4) its orbit. 



2. Observations and Data Reduction 



Subaru/IRCS Data - We imaged Fomalhaut on 2009 Aug 16 in the Mauna Kea 
J band filter behind the 0'.'8 coronographic spot with the Infrared Camera and Spectro- 
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graph (IRCS, iTokunaga et al.l Il998l ) and the AO- 188 adaptive optics system (full-width 
half- maximu m (FWHM) = (y/06; 20.53 mas pixel"^) in observing blocks bracke ting HR 8799 



obser vations ( ICurrie et al.ll2011al ) in angular differential imaging (ADI) mode ( iMarois et al. 
20061 ). Between the two observing blocks, we rotated the instrument by ~ 70° to keep the 
region covering Foma lhaut b visib l e on the dete ctor (5PA ~ 94°). Basic image processing 
follows steps listed in ICurrie et al.l (120101 l2011al Jb[) . 



HST/ACS Data and Basic Processing - We downloaded HST/ACS data for Fo- 
malhaut taken in the F435W, F606W, and F814W filters taken in 2004 and 2006 (Program 
IDs 10390 and 10598) from the MAST archive and processed with the OPUS pipeline. 

The distortion-corrected ACS data are pocketed by cosmic-ray hits. We flag most of 
these cosmic rays in individual subimages by identifying lO-o" outliers in concentric annuli 
roughly centered on the star in a median-filtered image with a moving-box length of 10 pixels 
We register images by subtracting the first image in the sequence by a 180° rotation of itself, 
identifying the centroid position that minimizes the residuals, and then cross-correlating 
this image with the others to derive relative offsets. We reform the subimages together into 
datacubes, mean-combining them with n-a outlier clipping to mitigate other cosmic-rays 
and hot pixels, iteratively choosing the exact settings to maximize the signal-to- noise of 
background objects (derived at the end of our image processing: next subject heading). 

PSF Subtraction - We explored a range of methods to subtract the stellar point- 
spread function (PSF), including reference PSF subtraction, roll subtraction /API-based PSF 
subtraction, LOCI ( "Locally Optimized Combin ation of Images" iLafreniere et al.ll2007l ). and 
adaptive LOCI or A-LOCI (ICurrie et al.l 120121 . T. Currie, 2012 in prep.). For the Subaru 
data, we obtain the best sensitivity using LOCI (A-LOCI yielded no improvement). For 
the F435W and F814W data, the small number of images or random noise sources reduce 
(A)LOCrs effectiveness. We obtain the smallest subtraction residuals by performing roll 
subtraction separately for each integration time and averaging the results. 



For the F606W data, A-LOCI provides the best results. The lCurrie et al.l version filters 
images by their degree of speckle correlation and determines the combination of algorithm 
parameters that maximize the SNR of point sources. Here, we add two components. First, 
we add a "moving pixel mask", where the algorithm is prevented from using, in the PSF 
construction, pixels lying within the annular region we want to subtract (the 'subtraction 
zone', see Marois et al. 2010b; Soummer et al. 2011). Second, we also include frames 
from the reference st ar Vega in our A-LOCI PSF construction (a "PSF reference libary"; 
Soummer et al.ll201ll ). 
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3. Results 

3.1. Detections, Astrometry and Photometry 

As shown in Figure [H we easily recover Fomalhaut b in all F606W data sets. With the 
noise defined in concentric annuli, the signal-to-noise ratios (SNR) are SNR = 16 for the 
3'.'0 2006 data (top-left panel), 9 for the I'.'S 2004 data (top-right panel), and ~ 8 for the 
2006 I'.'S data (not shown); the corresponding A-LOCl parameters are 5 = [1,1,1], Na = 
[125, 150, 125], g =[1,1,1], dr = [10,5,5], and Vcorr = [0.9,0.5,0.7]. Defining the significance 
of our detection within a 5-10 pixel annulus surrounding Fomalhaut b, we get SNR ~ 10.2, 
14.1, and 16.4 for the 2006 I'.'S data, S'.'O data, and the combined data (bottom panel) and 
SNR ~ 12.7 for the 2004 data. T hese estimates are about a factor of 2 higher than the best 



SNR estimates from iKalas et al.l (SNR ~ S; P. Kalas, pvt. comm.). The images include a 
background star (r, PA (2006) = 14'.'2, 203°) and two background galaxies (r, PA (2006) = 
lO'.'S, 135°; r, PA (2006) = 1S'.'7, 234°). We do n ot conclusively iden tify any other candidate 



companions at 5-cr, including at the location of iJanson et al.l (l2012[ )'s candidate nor within 
the debris ring. 



We also recover the FS14W detection reported by IKalas et al.l and report a new F435W 
detection of Fomalhaut b (SNR ~ 6) (Figure |2]). However, we do not detect Fomalhaut b in 
the Subaru/J band data (Figure |2]). 

To measure Fomalhaut b's brightness in the HST data and to quantify Subaru up- 
per limits, we u s e ape rture photometry and determine/correct for throughput losses as in 



Lafreniere et al.l ( 120071 ). Fomalhaut b is nearly unattenuated in A-LOCI reductions, a fact 
confirmed by comparing the F606W brightness of Fomalhaut b with photometry using ref- 
erence PSF subtraction. For the 2006 F606W data, we derive a brightness of m = 24.97 ± 
0.09. In the 2004 F606W data, we derive m = 24.92 ± 0.10. Fomalhaut b has not measur- 
ably varied in brightness between the two epochs, a result that appears insensitive to our 
choice of reduction method. Photometry at F435W and FS14W imply that Fomalhaut b is 
near zero color {'mF435w = 25.22 ± O.IS, nipsuw = 24.91 ± 0.20). Our 5-cr Subaru/J band 
limits are m > 22.22. 

We determine Fomalhaut b's position from the highest SNR data at each epoch (F606W), 
fitting a Gaussian profile (FWHM = 2.7 pixels). The astrometric uncertainties consider the 
intrinsic SNR, our precision in deriving absolute astrometric calibration, and differences in 
positions derived from the full range of PSF subtraction methods (reference PSF subtrac- 
tion, LOCI, A-LOCI). Fomalhaut b's position changes in the 2006 (2004) data by a[E,N] ~ 
[10, 10] mas ([15, 5] mas) depending on exactly how we do the PSF subtraction. The SNR 
for Fomalhaut b and background objects is insensitive (5SNR < 5-10%) to centroid offsets 
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< 13 mas (0.5 pixels) from the position that minimizes the subtraction residuals (see §2, 
paragraph 3). 

Our analyses yield a 2006 position of [E N]^^ = [-8.61 5, 9.352]" ± [0.016, 0.016]" and 2004 
position of [-8.598, 9.209]" ± [0.020, 0.014]. iKalas et al.r s astrometry shows good agreement 
with ours in the 2004 North coordinate but otherwise differs by ~ 0.5-1 pixels. Their 
astrometry also implies a larger space motion of ~ 0"18 vs. our ^ 0"14. 



3.2. Photometric and Astrometric Modeling 



3.2.1. SED Modeling: Planet Mass Upper Limits 



T o interpret the infrared non-dete ctions, we use t he pla net atmosphere models in lSpiegel and Burrows 
(120121 ) for 1-10 Mj companions and iBaraffe et al.l (120031 ) for 0.5-1 Mj companions. Our 



comparisons help to determine whether any planet masses are consistent with the combined 
upper limits and with the range of masses needed for Fomalhaut to sculpt the debris ring. 
Here, we con sider an age rang e of 100-500 Myr, encompassing age estimates from individual 
indicators in iMamajekl (120121 ). where the best-estimated age is ~ 450 Myr. 



Figure [3] compares several classes of models to the optical detections and the IR upper 
limits. The optical SED matches the spectrum of light from Fomalhaut (blue-green line) 
scattered by small dust. Optical photospheric predictions for the planet models are more 
than 5 magnitudes fainter than the 0.4-0.6 fim emissi on. Although the o ptical detections 
lie well above model predictions for gas giant planets, IJanson et al.l ( 120121 ) 's the 5-a IRAC 
upper limit restricts Fomalhaut b's mass to less than 2 Mj for even the oldest ages. The 
J-band data also restrict Fomalhaut b's mass to be < 2 Mj for younger ages, though useful 
limits at 400-500 Myr are entirely due to the IRAC data. 

However, the ground-based/Spitzer IR limits are still bright enough to be consistent 
with Fomalhaut b being a gas giant planet that sculpts the debris ring. The 5-cr IRAC 
upper li mits are comparable to the pr edict ed brightness of a 40 Myr-old, 1 Mj planet from 



both the ISpiegel and Burrows! (120121 ) and iBaraffe et al.l (120031 ) models. The IJanson et al. 
(I2OI2I ) upper limits are still ~ 3 (30) times brighter than a 0.5 Mj planet at (120) 500 Myr 



which is massive enoug h to sculpt the debri s ring: in fact, this low mass is favored by detailed 
dynamical simulations ( IChiang et al.ll2009[ ). 
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3.2.2. Nature of the Fomalhaut b Optical Emission 



Following iKalas et al.l (|2008[ ). a dust cloud producing Fomalhaut b's optical emission 



[mFQQQw,avg. = 24.95) at 119 AU requires a minimum cross-sectional area of cTp ~ l.SxlO^^/Qs 



1.3x10 



18 



m 



For spherically distributed dust, this implies a dust cloud radius of 



> 



0.0043 AU. Because Fomalhaut b is unresolved, the projected emitting radius of the cloud at 
F606W must be a < 0.5 FWHM or f« 0.26 AU. The timescale for orbital shear to broaden a 
dust cloud of initial ra.dius V r. into a region of radius a with orbital period P is ~ {tt/rc)P 
( iKenyon and Bromleyll2005l ). 



From this Keplerian shear alone, an unbound dust cloud at the separation of Fomalhaut 
b would then smooth out into a resolvable (a > 0.26 AU) clump in 56,000 yr and a circumstel- 
lar ring in ^ 25 Myr l/18th of Fomalhaut's age). Other dust grain removal mechanisms 
(radiation pressure, Poynting- Robertson drag) operate on similar/shorter timescales. Un- 
less we are fortuitously (and implausibly) identifying a very recent collision 
in a region with > 10 x lower optica l depth (and thus far loAv er collision frequencies) than 



the debris ring (cf iKalas et al.l 12005 
unbound dust cloud. 



20081 : iBoley et al.l 120121 ) , Fomalhaut b is not an 



3.2.3. Constraints on the Orbit of Fomalhaut b Derived from ACS Data 



Comparing the 2004 and 2006 ACS astro metry implies a deproject ed space velocity of 



3.7 ± 1.4 km s ^ or 33% slower than that derived by iKalas et al.l . If Fomalhaut b is 



a planet sculpting the debris disk, it s hould be nearly aps idally al igned with t he disk and 
have a space velocity of ~ 3.9 km s~^ (IChiang et al.ll2009l ). While iKalas et al.l 's velocity is 
~ 2-cr deviant from apsidal alignment {5.5~^qj km s~^), ours is consistent with that of an 
apsidally-aligned object within errors. 



Next, we compare Fomalhaut b's motion with the debris ring geometry, calculating the 
minimum distance between the companion and the disk midpoint (in angular separation). We 
explored a range of PSF subtraction methods for F606W 2004-2006 data sets to identify the 
combination that yields a high SNR image of the disk with minimal levels of self-subtraction. 
Based on these criteria, we use the a PSF subtracted image of the 2006 F606W data obtained 
under the 3'.'0 spot with Veg the reference star. 

To model the disk and compute the planet-disk separation, we explored two methods 
to assess the robustness of our results. Method 1 - We first select a wide annulus enclosing 
the disk emission, mask pixels with intensities roughly less than the 1-a noise floor (0.005 
cts/s) and determine a least-squares ellipse fit to the data. We perform these steps for 
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the raw image and a 10-pixel gaussian smoothed image, average the results, and report the 
differences in the results as the fit uncertainties. Method 2 - We divided the disk into 0.5° 
annular sections and fit the disk center at each section four different ways: 



• 1 A high-order polynomial fit to the surface brightness (SB) as a function of distance, 
taking the maximum value as the disk spine. 

• 2 A gaussian fit to the SB with the maximum of the gaussian as the disk spine. 

• 3 The peak pixel value 

• 4 The mean distance from all points > 90% the value of the peak pixel as identifying 
the disk spine. 



We then averaged the separately-determined disk radii, calculated x,y positions of the disk 
spine and uncertainties, and derived a least-squares ellipse fit. As with Method 1, we per- 
formed these steps on the raw image (5-pixel binning) and a gaussian smoothed image 
(3-pixel binning). 



Figure |3] (right panel) depicts our ellipse fitting from Method 1. Compared to lKalas et al. 



( I2OO5I ). we obtain similar values for the disk properties. From Method 1, we obtain 137.6 ± 
0.5 AU, 61.6 ± 0.1 AU, and 156.7 ± 0.2° for the disk major axis, minor axis, and position 
angle. From Method 2, we obtain 140.9 ± 1.7 AU, 58.0 ± 0.3 AU, and 156.6 ± 0.2°. From 
Method 1, Fomalhaut b appears to move closer to the disk by 0'.'026 (1'.'715 to 1'.'689), whereas 
Method 2 yields an even smaller deviation of 0'.'002 (1'.'712 to 1'.'710). This change is less than 
or equal to our quadrature-added astrometric errors for either the coordinate (cr2006+2004 = 
[0.026, 0.021]) and is not statistically significant. Thus, ACS astrometry is consistent with 
Fomalhaut b being in a ring- nested orbit and near apsidal alignment. However, uncertainties 
in the astrometry and disk fitting and the small number of data points do not yet rule out 
a ring-crossing orbit 0. 



^ KalasI ( 20101 ) reported preliminary results indicating that Fomalhaut b is at an implausibly wide separa- 
tion in the 2010 STIS data given their ACS astrometry. Ijanson et al.l (j2012l ) interpret this as evidence that 
Fomalhaut b is in a ring-crossing orbit. While a detailed analysis of the Fomalhaut STIS data is beyond 
the scope of this paper, our preliminary examination does not cast clear doubt on our results. We recover 
Fomalhaut b {msTiSfioew ~ 24.8 ± 0.2). Fomalhaut b's position is significantly less discrepant compared 
to the disk (« 0'.'1-0'.'15 ± 0'.'l)ifwe adopt the ACS ellipse parameters. Any uncorrected image distortion, 
clearly flagged as a problem by Kalasi will bias image registration by yielding erroneous determinations for 
where the diffraction spikes intersect (e.g. the stellar centroid). STIS astrometry of Fomalhaut b obtained 
at multiple epochs will better clarify its orbit (Kalas et al. 2012, in prep.). 
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Discussion 



Kalas et al.l ( 120081 ) argue that scattered starlight in a circumplanetary ring system con- 
tributes significantly to the F606W emission and planet thermal radiation dominates the 
F814W emission. They find the 0.6 iJ,m emission variable by .7-0.8 mag, e x plaina ble by 
accretion-driven emission, not dust emission. In contrast, iJanson et al.l ( 120121 ) argue 
that planet thermal emission cannot measurably contribute to the 0.6-0.8 /im emission else 
Fomalhaut b (if a planet) should have been detected at 4.5 fim. Based on the IR nondetec- 
tions, the large implied space velocity from the ACS data, the F606W varia bility, and the 
difficulty in a circumplanetary ring scattering enough starlight, iJanson et al.l instead argue 
that Fomalhaut b is a (semi-)transient dust cloud physically unassociated with any planet. 



We recover iKalas et al.r s detections at F606W and F814W and add a F435W detection, 
confirming that Fomalhaut b is a real object orbiting the primary. We cannot confirm that 
Fomalhaut b's F606W emission is variable at a level larger than our combined photometric 
errors (~ 0.15 mag), let alone by 0.7-0.8 mags, and thus do not find any evidence for 
accretion-driven emission nor for a semi-transient dust cloud. 

Based on the IRAC and (to a lesser extent) J-band upper limits, we agree that Fomal- 
haut b's optical emission contains no measurable contribution from a planetary photosphere. 
Thus, in the strictest sense of the w ord, Fomalhau t b is not yet a directly imaged planet. How- 
ever, a 0.5 Mj planet - favored by lChiang et al.r s simulations - is too faint to be detectable 
in any existing IR data set. Thus, the IR non- detections provide no compelling 
evidence against Fomalhaut b being the planet sculpting Fomalhaut's debris 
ring. 

Fomalhaut b's ACS positions appear to be consistent with those expected for a ring- 
nested orbit. Moreover, Fomalhaut b's implied space motion is significantly slower and much 
more consistent with that expe cted for an apsidally aligned object. Our astrometry then 
weakens the other argument in iJanson et al.l (120121 ) 's against Fomalhaut b being a planet 
sculpting the debris ring. 

Finally, Fomalhaut b is unlikely to be an unbound dust cloud since such a cloud would 
shear out into a resolvable clump on ver y short timescales and eventually into a circumstellar 
ring (see also iKuchner and Starkll2010l ). Alternatively, a dust cloud gravitationally bound 
to a planet can in principle explain the optical SED and persist for much longer times cales. 
Though dust confined to an optically-thick cir cumplanetary ring (IKalas et al.l l2008l ) may 
not scatter enough starlight (IJanson et al.ll2012l), a spherical cloud of circu mplanetary dust 
originating from ongoing planetesimal collisions (IKennedy and Wyattll201ll ) is still plausible. 



Fomalhaut b is a rare, possibly unique object in the primary's circumstellar environment. 
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Although the emission identifying it originates from dust, not a planet atmosphere, it just 
happens to be where a massive planet should be in order to sculpt the debris ring. If a planet, 
i t is (as yet) an indirect d etection, but this is a feature shared by the protoplanet LkCa 15 b 
( iKraus and Ireland! 1201 2| ) and radial- velocity/transit /mi crolensing-detected planets. Thus, 
Fomalhaut b is very plausibly (even 'likely') "a planet identified from direct imaging" even if 
images of it, strictly speaking, do not show thermal emission from a directly imaged planet. 



We thank Ray Jayawardhana, Paul Kalas, Mark Wyatt, Christian Thalmann, David 
Golimowski, Glenn Schneider and the anonymous referee for valuable comments. The STScI 
help desk staff graciously answered detailed queries regarding the HST/ACS data products 
and flux calibration. 
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Fig. 1. — A-LOCI processed images of Fomalhaut from 2006 (top- left; 3'.'0 coronagraph) 
and 2004 (top-right; I'.'S coronagraph) F606W data. The insets in the upper lefthand corner 
are centered on Fomalhaut b and have a box length of ~ 0'.'75 by O'.'S. Image combining the 
I'.'S and 3'.'0 2006 data rotated by 66° and enlarged with the inner r < 1'.'75 region masked. 



Fig. 2. — Roll- subtracted images of Fomalhaut in the F814W (top-left) and F435W (top- 
right) fihers from 2006 HST data and LOCI-processed J band/1.25 fim from 2009 Sub- 
aru/IRCS data (bottom). The image scale for the F814 and F435 data is the same as in 
Figure 1; the IRCS image focuses on a small, ~ I'.'O by 0'.'75 region centered on the measured 
and predicted positions for Fomalhaut b, with the 2009 (2004, 2006) circle size equal to the 
image FWHM (HST astrometric errors). 
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Table 1. Fomalhaut b Photometry 



UT Date Filter Processing MethodPhotometry (Vega) 



20090816 J 1 


> 22.22 (5-a) 


20060714-20F435W3 


25.22 ± 0.18 


F606W2,4 


24.97 ± 0.09 


F814W3 


24.91 ± 0.20 


20040926 F606W4 


24.92 ± 0.10 



Note. — Processing Methods: 1) LOCI, 2) Reference PSF Subtraction, 3) Roll Subtrac- 
tion, 4) and A-LOCl. 




Fig, 



(Left) Fomalhaut b photometry compared 



model spec tra from 



Spiegel and Burrowsl ( 20121 ) and photometry predictions from lBaraffe et al.l ( 20031 ) for 0.5-1 



Mj planets at 120-500 Myr and fit with a scaled mid-A star photosphere (blue-green spec- 
trum). The thick black lines and dotted black lines depict the transmission functions for the 
ground-based and Spitzer passbands, respectively. (Right) Ellipse fit to the debris ring from 
reference PSF-subtracted 3'.'0 F606W Fomalhaut data from Method 1 (units of counts/s). 
The yellow (blue) dot identifies the star (disk center). Like Figure 1 (bottom panel), the 
image is rotated by 66°. 



